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Abstract. We study the cosmological evolution of baryons in universes with and without 
primordial non-Gaussianities via N-body/hydrodynamical simulations, including gas cooling, 
star formation, stellar evolution, chemical enrichment from both population III and population 
II regimes, and feedback effects. We find that large /nl values for non-Gaussianities can alter 
the gas probability distribution functions, the metal pollution history, the halo baryon, gas 
and stellar fractions, mostly at early times. More precisely: (i) non-Gaussianities lead to an 
earlier evolution of primordial gas, structures, and star formation; (ii) metal enrichment starts 
earlier (with respect to the Gaussian scenario) in non-Gaussian models with larger /nl; (iii) 
gas fractions within the haloes are not significantly affected by the different values of /nl, 
with deviations of ^ 1% — 10%; (iv) the stellar fraction is quite sensitive to non-Gaussianities 
at early times, with discrepancies reaching up to a factor of ^ 10 at very high z, and rapidly 
converging at low z; (v) the trends at low redshift are independent from /nl: they are mostly 
led by the ongoing baryonic evolution and by the feedback mechanisms, which determine a 
^ 25% — 30% discrepancy in the baryon fraction of galaxy groups/clusters with respect to 
the cosmic values; (vi) non-Gaussianity impacts on the cluster X-ray emission or on the SZ 
effect(s) are expected to be not very large and dominated by feedback mechanisms, whereas 
some effects on the 21 -cm emission can be expected at early times; (vii) in order to address 
non-Gaussianities in the cosmological structure contest, high-redshift (z ^ 10) investigations 
are required: first stars, galaxies, quasars, and GRBs may be potential cosmological probes of 
non-Gaussianities. 



PACS numbers: 98.80 



1. Introduction 

Any cosmological model has to face two main concerns: the first is to explain the space-time 
evolution of the Universe, and the second is to explain the formation and growth of the cosmic 
"space-time content", i.e. the cosmic structures. 

Nowadays, it is possible to give good constraints on the first problem (Komatsu et al., 2010; 
Komatsu et al., 2011, and references therein): in particular, we know that the cosmological 
expansion is well described by the so-called ACDM model, in which the space-time is flat 
and its expansion rate at the present is Hq = 70km/s/Mpc (h = 0.7, if Hq is normalized 
to 100 km/s/Mpc). The density parameters describing the composition of the Universe are 
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(Komatsu et al., 2011) i^o.m = 0.272 for matter, and f2o,A = 0.728 for the cosmological 
constant A, often attributed to some kind of still unknown "dark-energy" (Maio et al., 2006, 
firstly investigated its effects on gas hydrodynamics and cosmic chemical evolution via 
detailed numerical simulations and found some impacts on the early gas cloud distributions). 
According to the cited observational data, the resulting total-density parameter is unity, 
^^o,tot = 1.0, but the only physically known fraction is very small, since baryons yields only 
l^o.b = 0.0456. 

The second problem is fundamentally related to the early "seeds" from which cosmological 
structures will develop. More precisely, clusters and galaxies are supposed to evolve from 
primordial matter fluctuations, that have a power-law primordial power spectrum, normalized 
via mass variance within 8 Mpc/ /i-radius sphere erg = 0.8, and with a spectral index n = 0.96. 
Then, in order to make statistical predictions for the structure growth, such as mass functions 
or bias (Wagner et al., 2010; D'Amico et al., 2011; Fedeli et al., 2011; Hamaus et al., 2010, for 
example), one has to assume a probability distribution of primordial fluctuations, around the 
mean density (p), as a function of the overdensity 6. According to the central-limit theorem, 
5 is usually assumed to be Gaussian. However, the validity of this assumption has often been 
questioned over the years and non-Gaussian models have been considered since long time 
(Peebles, 1983; Grinstein and Wise, 1986). Recent observational data do not rule out non- 
Gaussianities either (Komatsu et al., 2005; Komatsu, 2010). 

The most common way to parameterize non-Gaussian distributions (Desjacques and Seljak, 
2010; Verde, 2010, and references therein) is by introducing a non-linearity parameter, /nl, 
according to the following An^a^z (Salopek and Bond, 1990): 



In the previous equation O, $ is the Bardeen gauge-invariant potential, $l is the "linear" 
Gaussian part and the second term on the right-hand side is the non-linear correction, whose 
contribution is regulated by the /nl parameter. Since at any point the potential $ depends 
only on the value of $l, non-Gaussianities of this form are also called "local". Observational 
constrains on primordial matter fluctuations suggest that the primordial distribution could 
present some departures from the Gaussian one, with /nl in the range ~ 10 — 100 (Maio and 
lannuzzi, 2011, and their Table 2 for a collection of constraints). The existence of a non-null 
/nl would reflect on the growth of cosmological structures at any time and redshift (z). This 
means variations of the dark-matter halo mass function, mostly at large masses and/or at high 
redshift (Desjacques and Seljak, 2010; Verde, 2010; Lo Verde and Smith, 2011; Desjacques 
et al., 201 1, and references therein). Non-Gaussian effects could also be visible in the baryon 
history of cosmic objects (Maio and lannuzzi, 2011) and probably affect their gaseous and 
stellar components (Maio and Khochfar, 2011). The uncertainties on the baryonic evolution 
in non-Gaussian scenarios are mainly linked to the feedback effects which could strongly 
interfere with the original non-Gaussian signal. Besides the first attempt by (Maio and 
lannuzzi, 2011), there are, at the moment, no other investigations in this field, and quite 
a large lack of knowledge of the hydrodynamical interplay of the cosmic gas with non- 
Gaussianities. The ignorance of the detailed gas properties and of feedback effects could be a 
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strong limitation (Tescari et al., 201 1) when accounting for, e.g., metal abundances, Sunyaev- 
Zeldovic (SZ) signal (Sunyaev and Zeldovich, 1969; Sunyaev and Zeldovich, 1970), or X-ray 
emission from galaxy clusters in non-Gaussian models. Additionally, it might also influence 
the gas emission at high redshift (like the 21 -cm signal). 

Therefore, in the present study we analyze the principal features of the gas behaviour, during 
its condensation process within dark-matter haloes and address the effect on the baryon 
fraction in the Gaussian and in different non-Gaussian scenarios. We will focus on the cosmic 
gas and its distribution at high and low redshift, down to ^ = 0, by stressing the consequences 
from star formation, and the impacts for metal spreading, gas fractions, stellar fractions and 
baryon evolution within the formed cosmological structures. 

This work is organized as follows: in Sect. [2l we summarize the main properties of the 
simulations and of the physical implementations; in Sect. [3l we show our results about star 
formation (Sect. 13.11) . gas probability distribution functions (Sect. 13.21) . metal enrichment 
(Sect. 13.31) . and halo baryonic properties (Sect. 13.41) ; we finally discuss our results and 
summarize our conclusions in Sect. IH 

2. Simulations 

We use the simulations presented by (Maio and lannuzzi, 201 1), who performed the first study 
of the effects of primordial non-Gaussianities on the cosmic star formation, via self-consistent 
N-body hydrodynamical chemistry simulations. Beyond gravity and hydrodynamics, the 
code includes radiative gas cooling at high (Sutherland and Dopita, 1993) and low (Maio 
et al., 2007) temperatures, a multi-phase model for star formation, UV background radiation, 
and wind feedback (Katz et al., 1996; Springel and Hernquist, 2003, for further details). It also 
follows chemical evolution, metal pollution from population III (popIII) and/or population 
II (popll) stellar generations, ruled by a critical metallicity threshold of Zcru = 10"'' 
(Tornatore et al., 2007; Maio et al., 2009; Maio, 2009; Maio et al., 2010; Maio et al., 2011). 
The main agents of gas cooling (Sutherland and Dopita, 1993; Maio et al., 2007, e.g.) are 
atomic resonant transitions of H and He excitations that can rapidly bring temperatures down 
to ~ IC^K, when molecules become the only relevant coolants in pristine environments. 
In metal-rich environments, atomic fine-structure transitions are additionally included (Maio 
et al., 2007) and can further lower them, mostly at high redshift, since the formation of a 
UV background (Haardt and Madau, 1996; Dave et al., 1999; Gilmore et al., 2009) at lower 
redshifts heats the IGM and hinders gas cooling. In dense environments, stochastic star 
formation is assumed, cold gas is gradually converted into stars (Kennicutt, 1998), outflows 
take place, and feedback effects inject entropy into the surrounding environment, leading to a 
self-regulated star forming regime (Springel and Hernquist, 2003, for further detail). 
The simulations considered here assume a concordance ACDM model with total density 
parameter f2o,tot = 1-0, matter density parameter Vlo.m = 0.3, cosmological density parameter 
^o,A = 0.7, baryon density parameter i7o,b = 0.04, expansion rate at the present Hq = 
70km/s/Mpc, power spectrum normalization via mass variance within 8 Mpc//i-radius 
sphere ag = 0.9, and spectral index n = 1. We consider the cases: /nl= (Gaussian), 100, 
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1000, in cosmological periodic boxes with 100 Mpc//i a side, and with a gas mass resolution 

of ~ 3 X IO^Mq/L 

The stellar initial mass function (IMF) is a key ingredient for stellar systems, because it 
determines the temporal evolution and the relative distribution of stars with different masses 
and different metal yields. In the simplest approximation, it can be parameterized as a power 
law of the mass converted into stars (Salpeter, 1955) 

0(M.) oc (2) 

and normalized via the relation 

0(M,)dM, = l (3) 

,inf 

over the mass range [M^ „p]. Since the primordial IMF is quite uncertain, we bracket 

the possible scenarios by assuming both a top-heavy popIII IMF (normalized over the mass 
range [100, 500] solar masses, and with metal yields for PISN), and a (Salpeter, 1955)-like 
IMF (normalized over the mass range [0.1, 100] solar masses, and with metal yields for AGB, 
SNII, SNIa stars), as opposite limiting cases (Maio and lannuzzi, 201 1). The slope is fixed to 
X = —2.35 (Salpeter, 1955). For the popII-I stellar populations, we always adopt a common 
Salpeter IMF. 

3. Results 

In the next sections we show the main results of our analysis, starting, for clarity, from the star 
formation rate densities in the simulated volumes (Sect. 13.11) and considering the implications 
for the basic baryon properties, like probability distribution functions (Sect. 13.21) . and metal 
enrichment (Sect. [33]) . Then we will analyze the consequences on the cosmological structures 
(Sect. l34l) . exploring the evolution of their gaseous and stellar components (Sect. [3AT] 13.4.21 
and [3A3]) . 

3.1. Star formation 

We already mentioned that the N-body/hydrodynamical/chemistry simulations considered 
here have been presented for the first time in (Maio and lannuzzi, 201 1), who discussed the 
main effects of non-Gaussianities on the star formation history in the Universe. In the present 
work we aim at a more detailed study of the cosmological evolution of cosmic structures, and 
in particular of their visible components, gas and stars, within Gaussian and non-Gaussian 
contexts. Therefore, we start from the star formation rate densities,which are re-proposed, for 
sake of completeness, in Fig. [H for both the top-heavy popIII IMF (left) and the Salpeter-like 
popin IMF (right). The star formation rate density, p^,, quantifies the mass converted into 
stars, M^, per unit time, t, and unit volume, V, and can be aspressed as (Kennicutt, 1998) 
dp. ^ dM^ _p_ 
dt dtdv'^tff' 

where p is the gas density and tjf ~ l/y/Gp is the free-fall time. Therefore,the larger the 
gas density the larger the star formation rate density. As clear from the plots in Fig. [TJ the 
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Figure 1. Star formation rate densities as function of redshift for /nl=0 (solid lines), 
/f^L=100 (dotted lines), and /nl=1000 (dashed lines). Individual contributions from 
popIII regions (bottom lines) are also shown. In the left case the assumed popIII 
IMF is top-heavy (TH) in the mass range [100,500] M©, with a slope of -1.35 (Maio 
et al., 2010, have showed that the detailed value of the slope has little effects), while in 
the right case the assumed popIII IMF is Salpeter like (SL) (Maio and lannuzzi, 201 1, 
for further details). 

dependences on /nl are quite weak and it is not easy to distinguish among different cases, 
as long as /nl< 100. Only an extremely high (and currently ruled-out) value of /nl=1000 
would show appreciable differences with respect to the Gaussian case (/nl=0). 
Moreover, the global trends show that the Gaussian models give a fair representation of the 
star formation history, also with respect to /nl=100 non-Gaussian corrections. The main 
difference is found for the /nl=1000 case, which predicts a significantly earlier onset of star 
formation of almost lO^yr, and, thus, a correspondingly earlier enrichment process from the 
first stars and galaxies (see following sections). 

The adopted primordial popIII IMF has consequences on the popIII star formation. It 
determines a popIII contribution to the total star formation rate varying between almost two 
orders of magnitudes. More in detail, in the top-heavy popIII IMF case, the popIII star 
formation rate densities reach at most ~ 10^^ Mo/yr/Mpc^, for any /nl, and the relative 
contributions drop down below ~ 10^ at low redshift. In the Salpeter-like popIII IMF case, 
instead, the star formation rate densities achieved are ~ 10^^ Mo/yr/Mpc^, with a popIII 
contribution of less than ~ 10^^, even at low redshift, independently from /nl- 
We stress that the behaviour of the total star formation rate density at low z is practically 
unaffected by deviations from Gaussianity. 
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Figure 2. Gas probability distribution function at different redshift, z, as indicated by 
the labels, for the boxes with 100 Mpc/h a side, and a top-heavy popIII IMF. 



3.2. Gas distribution functions 

Star formation is the main agent leading the whole evolution of cosmic gas. Indeed, at 
very early times, when the Universe is younger than a few 10^ yr, there are no relevant 
star formation episodes, yet, and the gas has a purely pristine chemical composition. First 
stars synthesize heavy elements in their cores and eventually spread them around, during the 
final stages of their lives, by exploding as SN or PISN. These explosive events lead metal 
pollution in the surrounding regions and inject large amounts of entropy into the medium. 
Obviously, this has severe impacts (feedback) on the following generations, since the gas 
will be in very different conditions, with respect to the primordial, pristine gas. In order to 
illustrate and summarize the behaviour of cosmic gas during the cosmological evolution and 
its modifications due to non-Gaussianities, we consider the probability distributions, P{n), as 
functions of the number density, n: 

P(n) = 'M. (5) 
with f fraction of gas elements in the n density bin, and with normalization 

+ 00 

P{n)dn = 1. (6) 

-oo 

We display the gas probability distribution functions in Fig. |2l at different redshifts, for the 
runs having f]siL=0, /nl=100, and /nl=1000. The assumed popIII IMF is top-heavy. From 
the plots, it is evident that at early times (e.g. at 2; ~ 23 — 13) the gas distributions are slightly 
affected by non-Gaussianities, the gas is mainly at the mean density, and the environments 
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showing some deviations from the Gaussian behaviour are only the ones with over-densities 
~ 10^ — 10^, for /nl=1000. This basically reflects the features of the initial conditions, 
that, in the non-Gaussian cases, are biased towards higher densities. Within ~ 1 Gyr, a 
decrement in the differences of the high-density tails appears, with discrepancies going from 
~ 1 order of magnitude at 2; ~ 13 down to a factor of ~ 2 — 3 at 2; ~ 7.5. Moreover, below 
z ~ 7.5, the ongoing star formation events and the related feedback effects efficiently push the 
material far away from the star forming sites: this causes a larger spread in the distributions 
of the underdense material that is fed by gas expelled by the dense star forming regions. The 
feedback mechanisms cause gas mixing among the cold and hot phases, and the interactions 
between the inflowing and outflowing material increase the turbulent state up to Reynolds 
number of the order of ~ 10^ — 10^° (Maio et al., 201 1). These phenomena completely erase 
any imprints of non-Gaussianities (already by 2 ~ 3), for any /nl- This is an important 
conclusion, since it clearly shows that low-2; baryon evolution is affected by gas and stellar 
feedback in a much more important way than by deviations from the Gaussian shape of the 
underlying matter distribution (see further discussion in Sect. 13.4.21) . even when considering 
very large /nl- Gas emission at high redshift (e.g. the Hydrogen hyperfine-structure transition 
at 21 cm from early structures) is not supposed to change strongly for /nl=0 and /nl=100. 
For /nl=1000, whereas, there would be an enhancement due to the slightly larger amounts of 
gas at higher densities, but only within the first Gyr. 

The Salpeter-like popIII IMF cases present similar physical behaviours, with distributions 
which change accordingly to the different /nl values. 

More differences will be found in the metallicity evolution, as we will discuss in the next 
section. 

We note though (see also discussion in Sect. 131) the differences at early times in both gas 
distributions and star formation rates already suggests that primordial "visible" structures like 
first stars, galaxies, quasars, and gamma-ray bursts (GRBs) could be affected by the initial 
(Gaussian or non-Gaussian) dark-matter distribution. 



3.3. Metals 

Metal spreading is one of the principal signatures of the evolution of baryonic matter, and it 
leads crucial changes to the formation process of cosmic objects. 

To show the effects of different /nl values on metal evolution, in Fig. [3] we present the metal 
filling factor computed for the runs with fj^i,=0, /nl=100, and /nl=1000. This quantity 
represents the fraction of cosmic volume filled by metals and is computed from the outputs of 
our simulations, at each redshift, according to the equation: 

_J2imz,i/pz,i J2imz,i/pz,i 
Jv = 7 ~ 77 , (7) 

with i and j integers running over the number of the enriched particles and of all the SPH 
particles, respectively, niz particle metal mass, pz metal-mass density, m total particle mass, 
p total-mass density, and V simulation box volume. On the left panel, we show results from 
the simulations with a top-heavy popIII IMF, while, on the right panel, we show results from 
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Figure 3. Metal filling factors as function of redshift in the 100 Mpc/h side boxes for 
/j^L=0 (solid lines), /nl=100 (dotted lines), and /nl=1000 (dashed lines). Individual 
contributions from popIII regions (0 < Z < Zcrit, bottom lines) are also shown. In the 
left panel the assumed popIII IMF is top-heavy (TH) in the mass range [100,500] Mq, 
with a slope of -1.35 (Maio et al., 2010, have showed that the detailed value of the 
slope has little effects), while in the right panel the assumed popIII IMF is Salpeter- 
like (SL). 

the simulations with a Salpeter-like popIII IMF. In the plots, we also distinguish by metallicity 
regime to study the effects on the different stellar populations. The most evident feature is 
the convergence at low redshift of all the models. The most relevant differences are at high 
redshift, in correspondence of the different onsets of the first star forming events. 
In the top-heavy popIII IMF case (left panel), the first enrichment patterns are found 
immediately after the first star formation events, due to the extremely short lifetimes of 
massive stars and PISN (~ 2 x 10^ yr). Indeed, in the runs with /nl=0 and /nl=100 
the onset is at 2; ~ 16, and in the run with /nl=1000 it is at z ~ 22 (see left panel of 
Fig. [T]). Correspondingly, first metals are spread in the Universe at 2; ~ 15 for /nl=0 and 
/nl=100, and at 2; ~ 21 for /nl=1000. Given the large metal yields (Maio et al., 2010, 
for further discussions), the popIII filling factors become soon sub-dominant of some orders 
of magnitudes, for all the cases. We stress that the differences between the Gaussian case 
(/nl=0) and the /nl=100 case are negligible at any time. At lower redshift, when effects of 
non-Gaussianities are less relevant, all the models predict a total volume metal filling factor 
of ~ 10^^, and a popIII volume filling factor of ~ 10^'', four orders of magnitude smaller. 
When a Salpeter-like popIII IMF is adopted (right panel), the main changes are found at high 
redshift and are mainly due to the longer stellar lifetimes. In fact, they induce a larger delay 
of the pollution episodes, since the first metals are spread at z ~ 14 for /nl=0 and /nl=100, 
and at 2; ~ 18 for /nl=1000, i.e., about ~ 5 x 10'' yr, after the onset of star formation 
(consistently with the typical lifetimes of the first SNII explosions). Also in these scenarios 
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differences between /nl=0 and /nl=100 are not very appreciable and the asymptotic trends, 
led mostly by the popII-I regime, remain unaltered. 



3.4. Baryon fraction in the cosmic haloes 

After the general, global quantities previously studied, we investigate if and how the situation 
changes in the individual haloes found in the simulations, by paying more attention to their 
gaseous and stellar content. The formed structures are identified by a friends-of-friends (FoF) 
algorithm, with a linking length of 20% the mean inter-particle separation. In this way, we 
can determine total masses Mtot and basic properties of all the cosmic objects and of their 
different dark-matter, gaseous and stellar components. 



3.4.1. Gas and stars in the haloes In Fig. H the joint probability distributions P(/fe, Mtot) 
and P{fs,Mtot) for the baryon fraction, fh (upper panels), and for the stellar fraction, fg 
(lower panels), respectivelly, as a function of the halo total mass, Mtot, at redshift z = 0, are 
displayed: 

P{h,Mtot)= ^ dM^df ' 

and 

Pih,Mtot)= ^ HM^df • 

In both these relations, Nh^sampie is the total number of haloes sampled by the simulations, 
while is, in equation ([8]), the number of haloes within a given total-mass bin, M^^^, and 
within a given baryon-fraction bin, 4, and, in equation Q, the number of haloes within 
a given total-mass bin, M^^^, and within a given stellar-fraction bin, f^. With the previous 
definitions, the distributions are correctly normalized to 1. 
For the three /nl cases, we consider the runs with top-heavy popIII IMF. 
In Fig. m both baryon and stellar fractions have similar behaviours and no big differences are 
found. The baryon fraction, almost approaches the cosmic value at larger masses (a more 
quantitative analysis is performed in the next section) and exhibits some scatter (because of 
feedback mechanisms) in smaller structures, with peak values around fb ~ 10^^. The stellar 
fraction, fs, has a much larger scatter around the peak values of fs ~ 10"^, mostly in haloes 
with masses lower than ~ IO^'^Mq/Ji, and achieves asymptotic values of ~ 10^^'^ in the 
larger objects. In all these trends, though, no crucial impacts from different /nl are visible. 
The large scatter of and the even larger scatter of fs on galactic scales are strictly related to 
feedback effects that can strongly influence low-mass structures, because of their shallower 
potential wells. In fact, depending on the environment, gas heated by SN explosions in small 
haloes can be efficiently lost, and, on the other hand, more quiescent objects can be enriched 
by the gas expelled from nearby structures. Stellar scatter is probably connected to wind 
feedback, since each SPH particle undergoing star formation is eventually converted into a 
conlisionless star-like particle with a typical velocity of ~ 500km/s. Such velocities allow 
conlisionless stellar particles to escape from small-halo potentials and to end up into different 
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haloes. 

As a comparison, in Fig. [51 we show the status at redshift z ~ 7.3. In this case too there are 
no significant changes among the /nl cases considered. Nevertheless, in the /nl=1000 case 
there is slightly more baryonic (gaseous and stellar) mass in the larger objects, at masses of 
~ lO^^Mo/Zi: this reflects the growth enhancement at high redshift in models with large 
/nl values. The baryon fraction, f^, peaks at ~ 10^^, quite close to the cosmic value 
fb,cosmic = ^b/^m — 0.1333, for our choice of parameters. Since the star formation history 
at 2; ~ 7.3 has only started from just ~ 0.5 Gyr, the stellar fraction peaks at about one order 
of magnitude below the values at 2; = 0, i.e. at ~ 10^'^ — 10^^^. These results are almost 
independent from /nl, also in such high-z regimes, and keep some non-Gaussian signatures 
only at large masses. 

As a consequence, different /nl might slightly affect the masses of the haloes (see also 
Sect. 13.4.21) . but not significantly the gaseous and stellar fractions. In all the cases, the spread 
at low masses (at less than ~ 10^^ — 10^^ Mq/Zi) is independent from /nl, since it is similarly 
large and it becomes smaller and smaller at masses corresponding to galaxy groups or clusters, 
with the same trend for any /nl- This is essentially due to the capability of feedback in 
heating or expelling the gas from lower potentials, and enriching the closeby satellites. In 
larger objects, the material is more easily trapped and minor losses take place. In support of 
this, we will see in the next section that substantial amounts of baryons are lost during the 
bulk of cosmic star formation episodes, when feedback mechanisms become very important. 

3.4.2. Evolution of gaseous and stellar components To better check evolutionary differences 
due to the assumed value for /nl, vve follow the gaseous and stellar components within the 
most massive haloes found in the simulations, where also most of the star formation activity 
takes place. 

In Fig. |6l we plot masses (upper panels) and fractions (lower panels) of the baryonic (left), 
gaseous (center), and stellar (right) components, for the runs with /nl=0, /nl=100, /nl=1000. 
The assumed popIII IMF is a top-heavy one, while the dot-dot-dot-dashed line in the bottom 
left panel is the cosmic baryon fraction expected from the simulations. The baryon mass 
is about a few times IO^'^Mq/Ji at 2; ~ 10, and ~ IO^'^Mq/Ji at 2; = 0, hence this fairly 
represents the evolutionary path of a galaxy-cluster-like object. Also in this case, the mass 
evolution suggests that the model with the larger /nl=1000 predicts larger masses and earlier 
growth in all the components, while the models with /nl=0 (Gaussian) and /nl=100 do not 
have significant deviations. Looking at the temporal sequence, one sees that the major effects 
of non-Gaussianities are at high redshift (larger than z ~ 3), where differences reach up to 
a factor of ~ 2 — 3 (~ 10 in the very early phase) for the stellar mass, as highlighted in 
the bottom panels. Gas masses, instead, are affected by just ~ 1% — 10%, at most. The 
baryon fraction evolution shows that fb is roughly constant at z > 5 and traces the cosmic 
baryon fraction, fb,cosmic, quite well at early times. However, there are evident deviations at 
z below ~ 5 that reach the largest discrepancy between 2; ~ 1 — 3. The reason for that is 
immediately explained by Fig.fT] Indeed, at those z the star formation rate has its peak, and 
feedback mechanisms become more and more effective in removing material from the star 
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Figure 4. Baryon fraction, (upper panels), and stellar fraction, fs (lower panels), for 
all the halo total masses, Mtot, at redshift z = 0, in the simulations with /nl=0 (left), 
/p^L=100 (center), and /nl=1000 (right), color-coded according to their probability 
distribution (in logarithmic scale), P{fb, Mtot), and P{fs, Mtot), respectivelly. These 
results refer to the top-heavy popIII IMF cases. 




Figure 5. Same as Fig. |4l for redshift z = 7.33. 
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Figure 6. Redshift evolution of the masses (upper) and fractions (lower) of the 
baryonic (left), gaseous (center), and stellar (right) components, in the most massive 
halo, for /nl=0 (solid lines), /nl=100 (dotted lines), and /nl=1000 (dashed lines). 
Results refer to the top-heavy popIII IMF cases. The dot-dot-dot-dashed line in 
the bottom left panel is the cosmic baryon fraction expected from the simulations, 

fb,cosmic ~ ^b/^m- 





Figure 7. Same as Fig. |6l for a Salpeter-like popIII IMF. 



forming sites, and in regulating the star formation activity (as visible e.g. from the gas and 
stellar fraction trends in the bottom-center and bottom-right panels). 

Thus, star formation, also triggered by merger episodes, and the related feedback mechanisms 
are expected to be important sources of discrepancy between the cosmic baryon fraction and 
the one observationally estimated in galaxy groups and clusters (Giodini et al., 2009, e.g.). 
According to our findings, gas evacuation reduces of about 25% — 30%, mostly during the 
cluster assembly history, at z > 1, when smaller objects can loose more easily the material 
heated by feedback. E.g., an object with a baryon mass of ~ 10^^ Mq/H at redshift z = (as 
the one in Fig. (6]) is found to have a mass a few times smaller at redshift ~ 1 and more than 
ten times smaller at redshift ^ ~ 3, when its gas is being lost and the corresponding baryon 
fraction decreases. At z < 1, the gas results more tightly bound to the growing potential and 
some additional merging material, together with the decreasing star formation rate, sustains 
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Figure 8. Mass-weighted temperature (upper) and gas density (lower) profiles for the 
most massive halo found in the runs with /nl=0 (solid lines), /nl=100 (dotted lines), 
and /nl=1000 (dashed Unes). 



h. 

The role of non-Gaussianities in these regimes has little relevance, since, due to feedback 
effects on hydrodynamics, the high-2; enhancements of fs (due to large /nl values) gradually 
disappear and are completely absent below z ~ 3. The final values are always around 
/i, ~ 0.1, and/, - 0.03 - 0.04. 

In Fig. Ul a similar analysis for the Salpeter-like popIII IMF is done, but, besides the 
aforementioned differences at very early times (z ~ 8 — 10), the trends are still very similar 
and the physical interpretation does not change. 

3.4.3. Profiles Finally, in Fig. [8] the radial (r) profiles at 2; = for the gas mass-weighted 
temperature (Tmw) and density (p) of the most massive object are displayed in the three cases 
/nl=0, /nl=100j and /nl=1000 (for a top-heavy poplll IMF). The x-axis is r/r2oo, with r2oo 
the virial radiu^of each halo, equal to 1.064 Mpc/h, 1.131 Mpc//i, 0.736 Mpdh, for /nl=0, 
/nl=100, /nl=1000, respectively. The gas distribution in the haloes reaches overdensities of 
the order of ~ 10^ - 10^ (the critical density at ^ = is about 277.4 Mg/kpc^/i^), so the 
baryons are essentially decoupled from the cosmological frame. The global behaviour in the 
three cases is quite similar with a low-temperature environment (Tmw ~ 10^ — 10^ K) in the 
high-density regions (representing the cold innermost core) and an asymptotic behaviour at 
larger radii (where T^w ~ 10"^ K). We notice that, apart from the denser core in the /nl=1000 
model, no other clear trends with /nl are visible, and that the noisy inner part is more 

X The virial radius r2oo is defined as the radius at which the matter density is 200 times the critical density. 
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likely dominated by stochastic star formation feedback, gas cooling, non-thermal motions, 
etc. (Sarazin, 1988; Borgani and Kravtsov, 2009; Fang et al., 2009; Lau et al., 2009; Biffi 
et al., 2011). As a consequence, it seems that these structures have lost memory of the 
primordial initial non-Gaussianities and that their thermodynamical properties are mainly led 
by hydrodynamics and feedback mechanisms. These are typical examples of how the effects 
of different /nl values can be strongly mitigated or even cancelled by the more powerful 
effects coming from gaseous and stellar evolution within the growing haloes. 
We also point out that the drops in the temperature profiles at distances of a few per cent r2oo 
could be presumably alleviated by additional heating processes not explicitly considered here, 
like thermal conduction, cosmic rays or AGN feedback, and that are still widely debated in 
literature. An extensive study of all these different phenomena is beyond the aims of this work, 
but regarding non-Gaussianities we can safely state that they do not represent a crucial caveat, 
since they would add more degrees of stochasticity and would not affect our conclusions on 
the implications from different /nl parameters. 

4. Summary and conclusions 

In this work, we have considered the first large-scale N-body/hydrodynamical numerical 
simulations, dealing with non-Gaussian initial conditions, and including gas cooling, star 
formation, stellar evolution, chemical enrichment (from both population III and population 
II regimes), and feedback effects (Maio and lannuzzi, 2011). We have studied in detail the 
basic properties and evolution of baryons in the different scenarios, by paying attention to the 
effects on star formation (Sect. [3TI) on gas evolution (Sect. [3^ . on metal pollution (Sect. [33l) . 
and on the growth down to low redshift of the cosmological structures and their tgaseous and 
stellar content (Sect.[34l). 

The star formation (Fig. [U in Gaussian (/nl=0) and non-Gaussian (/nl=100, /nl=1000) 
models presents different behaviours, mostly at high redshift and for larger /nl, and, 
consistently with that, also the probability distribution functions (Fig.[2l) show more overdense 
gas in the /nl=1000 case and less in the /nl=100 and /nl=0 cases, to eventually converge at 
low redshift. 

The metal enrichment process is consequently affected (Fig. |3]), since metal pollution from 
dying stars starts earlier for /nl=1000 and later for /nl=100 and /nl=0- Different choices for 
the popIII IMF have effects at early times, as metals are very sensitive to the typical lifetimes 
of the first stars, but do not alter the trends for different /nl- 

Baryon, gas and star fractions within cosmological structures are very little affected by non- 
Gaussianities at low redshift (Fig.©, but might be quite affected at high redshift (Fig.[5l). In all 
the cases, typical values have quite a large scatter at low masses (< 10^^ M^/h at z = 0) and 
fb achieves almost the cosmic value at larger masses. More precisely, fb stabilizes a bit below 
the cosmic fraction with a weakly increasing trend for galaxy groups and clusters, while fg has 
a scatter of a few orders of magnitude below ~ 10^^ Mq/H and an almost constant or slightly 
decreasing trend at larger masses, where /^ ~ 10^^-^ (Giodini et al., 2009; Finoguenov 
et al., 2010, for observational studies at z < 1, and references therein). 
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Some clear differences due to the assumed /nl are more easily evident when following 
individual objects. In the case of the most massive clusters, gas mass evolution are affected 
by only ~ 1% — 10% at most, and stellar masses by up to a factor of ~ 2 — 3 at 2; > 4, both 
in the top-heavy popIII IMF scenario (Fig. (6]) and in the Salpeter-like scenario (Fig.|7]). These 
differences survive only at high redshift, though, since, the feedback mechanisms from gas 
and stellar evolution dominate the baryon behaviour at later times (Fig. [8]). 
The baryon fraction is significantly correlated to star formation and feedback, since we find 
discrepancies from the cosmic values (of ~ 25% — 30%) only when star formation becomes 
relevant (z < 4), independently from the primordial dark-matter distribution. 
Thus, it is likely that feedback mechanisms are the main responsible for washing out non- 
Gaussian signatures on the baryonic matter, by heating and mixing gas and by enhancing its 
turbulent state (Borgani et al., 2008). This would also mean that other gas properties, like 
temperature, density, velocity, can not keep vivid memory of the primordial /nl at lower 
redshift. The 21-cm emission from high-2 galaxies might still preserve some information, 
during the first Gyr, but the expected impact on, e.g., the X-ray emission from galaxy clusters 
or on the signal from the SZ effect(s) (Kitayama and Suto, 1997; Borgani et al., 1999; Sartoris 
et al., 2010) should be as small as the ones on the baryons, since they involve quantities, like 
electron fraction and temperture, that are strongly dominated by the feedback processes. 
As already mentioned in the previous sections, what emerges from our analysis is the need for 
high-redshift research, in order to investigate purely non-Gaussian effects. Moreover, thanks 
to available high-energy satellites (as e.g. Swi/q^l, Fermljjl Integra^, it is becoming 
possible to explore the 2; ~ 10 Universe also observationally (Cucchiara et al., 2011), 
and to impose constraints (based on observed data) on the GRB rate at high z (Campisi 
et al., 2011). Since the latter is strictly linked to the star formation history of the Universe 
and to its underlying cosmological model, GRBs might be suitable cosmological probes of 
non-Gaussianities, independent from cosmic microwave background (CMB) determinations. 
Also high-2; galaxies and quasars, could represent interesting probes of non-Gaussianities. 
The former are observed up to ^ ~ 10, and the latter are detected up to 2; > 6 together 
with many young, star-forming galaxies (Fan et al., 2001; Yan and Windhorst, 2004; Cool 
et al., 2006; Bouwens et al., 2007; Utsumi et al., 2010). Since these early structures are 
supposed to be tracers of high-density regions at high z (Djorgovski et al., 2003; Overzier 
et al., 2009), it is likely that their stellar content could be quite sensitive to non-Gaussianities, 
too. 

We conclude our discussion by summarizing as follows: 

(i) non-Gaussianities lead to an earlier evolution of primordial gas, structures, and star 
formation; 

(ii) the consequent metal enrichment starts earlier (with respect to the Gaussian scenario) in 
non-Gaussian models with larger /nl; 

(iii) gas fractions within the haloes are not significantly affected by the different values of /nl, 

§ http://heasarc.nasa.gov/docs/swift/swiftsc.html 

II http://fermi.gsfc.nasa.gov/ 

% http://science 1 .nasa.gov/missions/integral/ 
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with deviations of ~ 1% — 10%; 

(iv) the stellar fraction is quite sensitive to non-Gaussianities at early times, with discrepancies 
reaching up to a factor of ~ 10 at very high z, and rapidly converging at low z\ 

(v) the trends at low redshift are independent from /nl: they are mostly led by the ongoing 
baryonic evolution and by the feedback mechanisms, which determine a discrepancy in the 
baryon fraction of galaxy groups/clusters of ~ 25% — 30% with respect to the cosmic values; 

(vi) the impacts on the cluster X-ray emission or on the SZ effect(s) are expected to be not very 
large and dominated by feedback mechanisms, whereas 21 -cm from high-redshift galaxies 
might preserve more information related to /nl; 

(vii) in order to explore non-Gaussianities and baryonic structure formation, high-redshift 
{z ~ 10) investigations are required, with first stars, galaxies, quasars, and GRBs being 
optimal fields of study, and potential cosmological probes of non-Gaussianites. 
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